Self-interfering matter-wave patterns generated by a moving laser obstacle in a 
two-dimensional Bose-Einstein condensate inside a power trap cut off by box potential 
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We report the observation of highly energetic self-interfering matter-wave (SIMW) patterns gen- 
erated by a moving obstacle in a two-dimensional Bose-Einstein condensate (BEC) inside a power 
trap cut off by hard-wall box potential boundaries. The obstacle initially excites circular dispersive 
waves radiating away from the center of the trap which are reflected from hard-wall box bound- 
aries at the edges of the trap. The resulting interference between outgoing waves from the center 
of the trap and reflected waves from the box boundaries institutes, to the best of our knowledge, 
unprecedented standing wave patterns. For this purpose we simulated the time dependent Gross- 
Pitaevskii equation using the split-step Crank-Nicolson method. The obstacle is modelled by a 
moving impenetrable Gaussian potential barrier. Various trapping geometries are considered. 



I. INTRODUCTION 

Recently, there has been considerable interest in the 
investigation of the effects of a moving obstacle (or ob- 
ject) through a Bose-Einstein condensed system in vari- 
ous trapping geometries [IKS]. The obstacle is a poten- 
tial barrier generated by a Gaussian laser beam [3J [5], 
which can be repulsive or attractive depending on the 
vifavelength of the beam. Previous work has also shown 
that the motion of this obstacle, whether linear or rota- 
tional, causes excitations leading to strikingly interesting 
phenomena such as vortices [H flMIi] , solitons [TU [TS] , 
crescent vortex solitons |16j , and dispersive waves accom- 
panying the solitons [T7]. In addition, there has been 
growing interest in the investigation of Cerenkov radia- 
tion [7j i8j and waves generated by an obstacle moving 
with supersonic velocity [H |6] . However, to the best of 
our knowledge, none of the previous literature considered 
the effects of adding a hard-wall to the magnetic trap on 
the dynamics of the BEG excited by a moving obstacle. 
It is well known, that vortices can be excited by moving 
obstacles inside a trapped BEG. However, the parameters 
used in the present investigation arc lower than what is 
required to obtain significant vortex features, such as the 
vortex-antivortex pairs observed by Jackson et al. [1 . 

One goal of this work is to check for a possible exis- 
tence of vortices in our systems below. Our motivation 
for this stems from the vast literature on vortices in BEGs 
confined by various traps, with interesting and exciting 
phenomena both in theoretical [TUl El [T5li24j and exper- 
imental domains p5H28j . 

One way to explore trapped BEGs excited by a mov- 
ing obstacle is by the time-dependent Gross-Pitaevskii 
equation (TDGPE) whose solutions are known to be 
solitons "29] depending on the interactions and certain 
boundary conditions; the nonlinearity in the GPE being 
the cause for the appearance of these excitations. The 



role of the nonlinearity in the generation of solitonic so- 
lutions for the nonlinear Schrodinger equation has been 
studied as early as 1983, when Nonnenmacher and Non- 
nenmacher |29j found, that for repulsive interactions no 
solitons are observed; unlike the case for attractive inter- 
actions. 

In this work, the split-step Grank-Nicolson method 
|30) is applied to numerically solve the two-dimensional 
(2D ) TDGPE, involving a harmonic or (power-law) PL 
trap plus a moving obstacle. The trap is cut off by hard- 
wall box-potential (HWBP) boundaries serving as reflec- 
tors of matter waves. In fact, the box potential, along 
with other hard-wall trapping geometries, has been used 
previously by Ruostekoski et al. [31 , who explored the 
interference of a BEG in a hard-wall trap plus the non- 
linear Talbot effect. We shall comment on their work in 
the discussion section. According to Ruostekoski et al., 
the hard-walls can be realized experimentally by a blue- 
detuned light sheet. Gonsequently, the trap is limited by 
the hard walls and an expanding BEG can be reflected 
off these walls producing a self-interfering matter-wave 
(SIMW) field. Essentially, then, we propose an experi- 
ment similar to that of Ruostekoski et al., in which we use 
a 2D^ square box potential. But instead of switching off 
the isotropic magnetic trap in order to allow the BEG to 
expand, we excite the BEG by a moving obstacle in order 
to generate dispersive waves (DWs) that travel towards 
the hard wall in order to get reflected. Thus, if no DWs 
are emitted, no SIMWs can be obtained. However, our 
2D square box is achieved numerically by a different 
method than that of the latter authors. The main mo- 
tivation of this work, then, is largely based on the latter 
work of Ruostekoski et al. |31j . Additional motivation 
stems from a previous investigation by Jackson et al. [T] , 
who simulated the motion of an "object" through a di- 
lute BEG in a 2D harmonic trap, where the object (or 
obstacle) corresponds to a laser-induced vortex potential. 



Our goal is to explore the same system of Jackson et ai, 
at a lower mean-field interaction strength Af but with 
an additional box-potential trap. We essentially take a 
closer look at the dynamics of the density surrounding 
the central BEG at distances larger than ~ 4 trap lengths 
from the center of the trap in the event of matter-wave 
reflections. Another motivation comes from the work of 
Hakim [17', who investigated the superflow of a nonlin- 
ear Schrodinger fluid past an obstacle. He found that 
the obstacle repeatedly emits gray solitons which prop- 
agate downstream; whereas DWs propagating upstream 
are emitted at the same time. Since the calculations of 
Ref. [T7\ were conducted in ID and outside a trap, the 
importance of our current investigation lies, then, in ex- 
ploring these systems in 2D and further inside various 
trapping geometries inside a hard wall potential. It is, in 
fact, the DWs which cause the phenomena obtained in 
our current study. To this end, we chiefly aim at present- 
ing a phenomenological investigation of SIMW patterns 
generated by a moving obstacle inside a 2D BEG in 
various trapping geometries surrounded by a HWBP. 

Several questions are tackled: First, can we observe 
DWs in 2D^ harmonically or nonharmonically-trapped 
BEGs excited by a moving obstacle, such as was observed 
by Hakim T7] for the case of a ID uniform Bose gas? 
Second, by adding the HWBP to the harmonic or PL 
trap, what new features are observed when the latter 
matter waves are reflected from the hard walls? Will 
they be similar to Ref. [31 j? Third, how does the energy 
of the system behave with time? Fourth, what about the 
momentum density dynamics? 

Further, additional simulations of the ID GPE are 
conducted for the same systems above in an attempt to 
reproduce the dispersive waves observed by Hakim llTj. 

The purpose is to compare these ID^ disturbances with 
the 2D^ disturbances in order to understand their nature 
in 2D. 

Our key results are as follows, i) It is found that in 

a 2D^ harmonic or PL trap, highly-energetic DWs can 
be emitted from the moving obstacle. These DWs are 
highly energetic, capable of surmounting the PL poten- 
tial barrier to far distances from the trap center. Upon 
reflection of their wave fronts from the HWBP, their in- 
terference with incoming DWs produces at some point 
perpendicular, 2D density-stripes in the matter-wave 
field, largely influenced by the artefacts of the trapping 
geometry. This phenomenon has not, to the best of our 
knowledge, been reported elsewhere to this date; ii) the 
accompanying energy dynamics of these systems display 
an oscillatory pattern indicative of a possible presence of 
solitons; iii) when the trapping strength of the PL is in- 
creased, the energy otherwise expended on the creation 
of the SIMWs is channeled to excite a larger number of 
solitons inside the trap, eventually leading to a "soliton 
gas"; iv) the condensate dynamics reveal further an os- 
cillatory behaviour signalling that particles are excited 
out of and deexcited back into the BEG. 



The paper is organized as follows. In Sec|nj we briefly 
describe the method used. In Sec |III[ we present our 
results, and in Sec |IV| we discuss them. The conclusion 
is presented in Sec|V| 



II. METHOD 
A. Time-dependent Gross-Pitaevskii Equation 

In this work, we apply the split-step Grank-Nicolson 
method [30] to numerically solve the 2D^ TDGPE given 
by 
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where N is the number of particles, g = AttTi as/m is 
the interaction parameter, with Os the s-wave scattering 
length in the low-energy and long-wavelength approxi- 
mation, m is the mass of the boson, and Ti is Planck's 
constant. For the present purposes, V{x,y]T) is taken 
as a potential which varies with time. According to Mu- 
ruganandam and Adhikari (MA) J30], the TDGPE can 
then be recast into a dimensionless form using units of 
the trap such that 
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where t = uj^oT is a unitless time, V{x,y]t) 
V{x,y;t)/{huJho) and x -^ x/t, y -^ yji with i 
^JhJJjnujho) the trap length. The parameter J\f is 



N = 



2Na, 



(3) 



A being a parameter describing the width of the ground- 
state wave function 0o(z) = (A/7r)-'^/'*exp(— A2;^/2) in 
the z-direction, which was then integrated out from the 
3D^^ wave function so as to get the 2D^^ form Eq.([2| 
above. That is, we consider the z-direction to have har- 
monic confinement of infinite strength. Note further 
that MA rescaled the wave function using Lp{x, y; t) = 
%/F*(r;T). 



B. Trapping potential and moving obstacle 



The time-dependent trapping potential is composed 
of two parts: a static external general 2D^ PL trap 



VpL{x,y), and a moving potential (MP) barrier (or ob- 
stacle) with time VMp{x,y;T) which is generated by 
sweeping a blue-detuned laser beam through the BEC 
(see, e.g., [52] )• Hence, V{x,jj;t) = VpL{x,y) + 
VMp{x,y;t). The PL potential VpLix,y) is taken to be 
of the form 



VpLix.y) = a{\xr +^\yn, 



(4) 



where Pi, P2 > can be any positive numbers, k is an 
anisotropy parameter, and a is the strength of the PL 
potential, taken here to be 1. Thus for pi = p2 = 2 we 
have the usual harmonic trap, whether isotropic (k = 1) 
or anisotropic (k > 1 or < 1). Both a and k are unitless, 
since x and y are unitless as well. The goal of using 
Eq. Q is to explore the effect of different curvatures of the 
external trapping potential on the dynamics of a trapped 
BEC when excited by a moving obstacle. 

Next, in order to excite the BEC, we follow Jackson 
et al. pi] and Bao and Du [33] and introduce the obstacle 
in the form of a velocity-dependent Gaussian potential 
given by 



VMp{x,y;t) = Aexp[~/3x^ - /3(y - vt)^ 



(5) 



A being the amplitude, v the velocity in the y— direction, 
and /3 the exponent determining the width. 

We further cut off the harmonic or PL trap by a 
HWBP, which can be numerically achieved by forcing 
the wave function of the system to vanish at the hard 
walls. This can be realized by imposing the following 
boundary conditions 



(p{x,~Ly/2-t) = ip{x,Ly/2;t) = 
ipi-L,/2,y;t) = ip{L,/2,y;t) = 0. 



(6) 



Thus the overall trapping geometry is then described by 
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We take L^ = Ly = 20 in trap units. Note that this is 
different from the way the box potential was constructed 
in Ref. [3T] , where a set of large-amplitude Gaussians was 
used. In essence, we mimic here an "exact" hard-wall 
potential where the BEC cannot tunnel —even slightly- 
through the wall. 



C. Energy functional 

The energy is evaluated by the well-known energy 
functional in trap units 
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Since the above functional includes the moving obstacle 
Eq.((5]), the energy is time-dependent. In fact, an oscilla- 
tory pattern is revealed in Sec III E 



D. Momentum density 

In an attempt to investigate the zero-momentum, as 
well as higher-momentum states, we set out to compute 
the momentum density of the system given by the Fourier 
transform 

PFT{kx,ky;t) = ^-^ / / pix,y;t)e'''''''+'''yydxdy, 

(9) 
where p{x,y;t) = |0(r;i)p and {kx.,ky) is a 2D mo- 
mentum state. The goal is to explore the dynamics of 
the condensate fraction as well as the momentum den- 
sity distribution. 

E. Numerics and initial conditions 

We applied a previously-written Fortran 77 code [SD] 
to numerically solve Eq.(2]) in real time, using the split- 
step Crank-Nicolson method. For our purposes, we mod- 
ified the code by including the PL potential Eq.Q and 
the moving Gaussian potential Eq.([5|, instead of the 
anisotropic harmonic trap already present in the code. 
The simulations were conducted on a grid of 400 x 400 
square pixels, the edge of each having the size of a step, 
0.05. The time step was chosen to be A = 0.001. The 
code has a part which initializes the system and a tran- 
sient part during which the system evolves. In the ini- 
tialization step, the nonlinearity is introduced gradually 
(adiabatically) by a stepwise introduction of the inter- 
action parameter A/" at a rate of dAf/dt = Af/{Nstp/^), 
where Ngtp is the number of time steps A during which 
the system is initialized. In essence, we modified the ini- 
tialization part to include a gradual introduction of the 
Gaussian potential Eq.(l5]) by a similar gradual increase 
of its parameters. In fact, we explore in this paper two 
different initialization conditions via which the obstacle 
is introduced into the system along with the nonlinear- 
ity. It will be later shown, that the generation of DWs 
depends strongly on the initialization conditions of the 
system. We find that one condition [(i) below] leads to 



the generation of DWs, whereas the other [(ii) below] 
does not. These conditions are as fohows: 



i) Simuhaneously with the nonhnearity-introduction, 
the system is initiahy subjected to a moving obsta- 
cle whose amplitude A and velocity v are gradu- 
ally increased with time from zero until they reach 
the desired specific values in this simulation. The 
corresponding rates are dA/dt — A/{NstpA) and 
dv/dt — ^/(TVgtp A), respectively. At the end of the 
initialization, the obstacle will have left the center 
of the BEC. Subsequently, a second equivalent ob- 
stacle is again abruptly switched on at the center 
of the BEC, with the same parameters reached at 
the end of the previous initialization. Then it is set 
into motion with speed v until it exits the trap. 

ii) Simultaneously with the nonhnearity-introduction, 
stationary obstacle is applied in the initialization 
process, whose amplitude is gradually increased at 
the rate A/{Nstp^)- After initialization, the obsta- 
cle is set into motion with speed v until it exits the 
trap. 



For the interactions and obstacle we used the parameter 
values N' = 5, A = 30, P = 3, and v = 2. Further, we 
used Nstp — 42000 such that the rates dA/dt and dv/dt, 
using A — 0.001, are 0.714 and 0.048, respectively. Sim- 
ilarly, Af is gradually increased at a rate of A/'/(A^stpA) — 
0.119. For ^^Rb atoms with some of the same parameters 
used in Ref. [31 , i.e., an s-wave scattering length a^ = 5.4 
nm and a trapping frequency uJho = 2tt x 25 Hz, the trap 
length iia -^ 2.16 x 10'* A. If we choose A = 100 so that 
the width of the ground state (^o(-2) becomes extremely 
small along the z— axis, one gets A'' ~ 40 particles. If we 
box length is taken 2L — 15/im, then the density of our 
system is n ^^ N/{4L'^) — 4 x 10^°m~^. Consequently, 
our systems are in the dilute regime na^ ^ 10^^. Ref.[3T] 
used particle numbers of the order ~ lO'^, and thus we 
aim at checking the effect of using a much lower particle 
number. 

Experimentally, then, i) one initially sweeps a very 
weak laser beam whose intensity and velocity are grad- 
ually increasing. After the initialization, when the ini- 
tial laser beam has left the BEC, a second laser beam is 
abruptly switched on at the center of the BEC. Then it 
is set into motion with a constant velocity or accelera- 
tion. The second laser beam has the same characteristics 
as its predecessor at the end of the later initialization, 
ii) One subjects the system to a stationary laser beam 
whose intensity is gradually increasing with time. After 
the intensity has reached a maximum, the laser beam is 
set into motion with a constant velocity or acceleration. 

We therefore propose experiments in which one could 
try to verify the influence of the latter conditions on the 
generation of DWs. 



III. RESULTS 

In what follows, we present the results of our simu- 
lations. Our major findings are as follows. The mov- 
ing obstacle inside the BEC generates, in addition to 
vortices [TJ [T0Vll4j . 2D^ circular DWs radiating away 
from the moving obstacle. As their wave fronts reach 
the hard walls of the box-potential, they are reflected 
backwards upon which they interfere with incoming cir- 
cular DWs. This interference creates interesting self- 
interacting matter-wave (SIMW) patterns. On increasing 
the powers pi and p2 of Eq. Q to values beyond 2.0, these 
SIMWs are suppressed. The suppression prevents the en- 
ergy loss from being carried away by the DWs and saves 
the energy for another purpose. In the present case, this 
energy is directed into another channel in which it is ex- 
pended on exciting a larger number of solitons inside the 
trap. Similarly, on increasing the anisotropy of the har- 
monic or PL trap, DWs travelling parallel to the minor 
axis of the trap, are also suppressed. We have reasons 
for speculating the presence of real solitons in our sys- 
tems as opposed to only solitonlike structures reported 
earlier in Ref . [31j . This is largely supported by the os- 
cillatory behaviour of the energy dynamics indicative of 
soliton generation |34j . By inspecting the density struc- 
ture at the center of the BEC at some time, the solitons 
turn out to be of a vortex structure. Finally, the mo- 
mentum distribution of the system reveals the presence 
of a zero-momentum condensate density and the conden- 
sate fraction oscillates with time. We would like to draw 
the attention of the reader, that all the upcoming re- 
sults have been recorded after an initialization process of 
length t — 4.2. Thus, when for example a time t = 3 
is mentioned, it refers to the transient time of the sim- 
ulation. The total simulation time is then 7.2 for this 
particular case. However, we do not record results in 
the initialization stage, as we would like to inspect post- 
initialization phenomena. The following results were ob- 
tained using the initial condition (i) in Sec II E It is only 
in Sec 



SecHI 



III F that we discuss the effect of condition (ii) in 



A. Harmonic trap with an obstacle of constant 
velocity 

Figure [ij displays a 2D^ density plot obtained at a 
transient time of t = 3 for A/" = 5.0, pi — p2 = 2 (har- 
monic trap), and n — 1.0. Fig. [2] displays the corre- 
sponding density map. For the obstacle we used A — 30, 
/3 = 3.0, and v — 2.0. The intriguing result is that four, 

mutually perpendicular, 2D SIMW stripes arise. Two 
of them propagate parallel to the x— axis, the others par- 
allel to the y— axis along the motion of the obstacle. In 
the area surrounding the latter cross-like pattern, some 
residual faint interference pattern is also present. We 
shall relate this phenomenon to an artefact of the trap- 
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FIG. 1: Density of a 2D^' trapped BEC excited by a moving 
obstacle, Eq. ([5]l, as obtained by a Crank-Nicolson simulation of 

the 2D^1 time-dependent GPE, Eq.H, at t = 3.0. The parameters 
used are A^ = 5.0, pi = P2 = 2.0, ana k = 1.0. As for the obstacle, 
Eq.(l5]l, we used A=30, 13=3.0, and d=2.0. The Fortran program 
usedwas that of Adhikari and Muruganandam |30) . Lengths are 

in units of the trap [a(,o = \/h/{mujfia) ]■ 




FIG. 2: Density map of Fig. IT] Lengths are in units of the trap 

[aho = \/^/imu}ho) ]■ 



ping geometry discussed later in Sec. |IVB| It is noted 
that the SIMW stripes generated in the regime y > are 
perturbed by the obstacle causing the density depression 
(small valley) in the upstream direction, behind the cen- 
tral BEC peak at t = 3. Further, upon a closer inspection 
of the area surrounding the small moving "valley" in the 
upper part of Fig. [2] a pair of "holes" can be seen iden- 
tifying a pair of vortices. This valley displays also a com- 
plex structure which could be due to vortex excitations. 
It should be recalled that the motion of the obstacle is 
in the positive y— direction. For further demonstration, 
see the movie in the supplementary material section [35] 
corresponding to the time evolution of \(p{x, y; t)p in Fig. 
Ill Upon a careful inspection of the movie, one can see a 



FIG. 3: Same as in Fig. IT] but for an anisotropic trap with k = 4.0 
[Eq.Ml]. Lengths are in units of the trap [ajj^ = i/h/{mujfia) ]■ 



solitonlike structure in the matter-wave field upon the re- 
flection of the BEC from the hard walls, similar to what 
has been reported by Ref. |31j. 

Figure |3| displays the 2D^ density for the same 
parameters as in Fig. [l] except that the trap is now 
anisotropic with k — 4, with x the major and y the 
minor axis. Because of a stronger confinement along 

the J/— axis, 2D^ SIMW stripes are absent in that di- 
rection, as outgoing matter waves from the obstacle are 
suppressed, unable to tunnel through the potential bar- 
rier. Rather, they are observed to emerge only along the 
positive x— axis. There is another matter wave emerging 
parallel to the 2D^^ SIMW in Fig. [s) which has a larger 
wavelength and therefore less energy than the accompa- 
nying SIMW. 

In Fig. |4J a magnified view of a section of the SIMW 
stripes appearing in Fig. IT] is presented. Apparently, the 
SIMWs eventually propagate almost without a significant 
change in width or amplitude. Further, they display a 
peculiar complex structure. Another magnified view of 
a section of the SIMWs of the anisotropic case of Fig. 
|3] is displayed in Fig. [5] This time, two separate wave 
structures are revealed. One is the SIMW centered along 
the y = axis; the other is a decaying matter wave of a 
larger wave length centered along the y = —4 axis. The 
latter second wave decays after a certain distance (^ 10 
trap lengths) from x — and is therefore less energetic 
than the DWs. 



B. Harmonic trap with an accelerated obstacle 

The purpose of this section is to demonstrate that 
2D SIMW stripes can still be generated by applying 
an accelerated obstacle (AO). For this purpose, a variant 
of the Gaussian potential [Eq.(l5])] is introduced, where 
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FIG. 4: Magnified section of the SIMW stripes appearing in 
Fig. IT] showing how these waves propagate. Lengths are in 
trap units [auo = ^/h/(muJ^] . 
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FIG. 6: Same as in Fig. IT] but with an accelerated obstacle with 
a = 1.0 [cf. Eq.lToll]. 



FIG. 5: Same as in Fig. |4] but for Fig. [Sjwith «:=4.0. There 
are now two separate wave structures. One centered along 
the y = Q axis which is a SIMW; the other along the y = —4 
axis. Lengths are in trap units [aho ~ W^/i'mujho)]- 



the velocity term is replaced by an acceleration term 

(l/2)ai2: 



Vao = ^exp 



-I3x^ - (3[y- -at 



(10) 



where a is the acceleration of the potential. Another goal 
is to display the effect of the "force" of an obstacle on 

the dynamics of a 2D trapped BEC. The obstacle has 
a "mass", as has been outlined earlier by Astrakharchik 
and Pitaevskii [5] , and might therefore be able to produce 
a shock wave when accelerated. 

Figure [6] displays the resulting density plot at i = 3 
for the same parameters as in Fig. IT] except that we are 
using the AO Eq.([lO| with a = 1 instead of Eq.([5]). It 
can be seen that, again similarly to Fig. IT] a perpendic- 
ular set of 2D^ SIMW stripes is observed. As a result, 
one concludes that an accelerated obstacle is also able to 
generate 20^^ SIMW stripes in a 2D^1 trapped BEC. 
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FIG. 7: Same as in Fig. [T] but for pi = p2=2.2 at t = 2.4. 
Lengths are in units of the trap [a^o = \/'hl(rrujJho)\- 



Nonharmonic potentials: effect of trap 
curvature 



In this section, we explore whether one can still obtain 
2D SIMW stripes if one uses a stronger PL trap with pi , 
P2 > 2. It was found that for a choice of pi — "pi = 2.2, 
one still obtains these stripes similar to Fig. [T] as dis- 
played in Fig. W\ However, a closer inspection of the 
latter figure reveals a pair of parallel stripes propagating 
downstream along the negative y— axis. The parameters 
are the same as those used in Fig. fTl except for pi and 
P2- The intensity of the 2D SIMW stripes is, however, 
smaller parallel to the a;— axis than the y— axis. On go- 
ing up to larger powers, pi = P2 > 3, the 2D SIMW 



W{x,y]t)\ 



AT = 5, yl = 30, K = 
^ = 3, i' = 2, t = 2 
Pi = P2 = 3 




X {aho) 



FIG. 8: Same as in Fig. [T] but for pi = p2 = 3.0 at t = 2.0. 
Lengths are in units of the trap [aho ~ \/T^/{mu]ho)\- 
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Pi = P2 = 4 




FIG. 9: Same as in Fig. [T] but for pi = p2 = 4.0 at t = 2.0. 
Lengths are in units of the trap [aho ~ yT^/i'mijJho)]- 



stripes are suppressed as displayed in Figs. [8]and|9] The 
curvature at the edges of the trap plays a crucial role in 
the suppression of the SIMWs as explained later in Sec. 
lYl 



D. One-dimensional simulations 

In this section, we present additional results of the 
Crank- Nicolson simulation of the ID^ GPE, i.e., the 
ID^l version of Eq.lB: 



1^ 
'2 5x2 



V{x;t) + A/'|¥>(a;;i)|^ 



'dt 



(p{x;t) = 0, 
(11) 



where 




X (aho) 

FIG. 10: One-dimensional SIMWs generated by a ID^^ 
Cra nk- Nicolson simulation of the time-dependent GPE 
Eq.([ll|. The system is a ID^ hard-sphere Bose gas of 
Af — 5. The vortex potential parameters are A = 10, 
13 — 3, and v = 2. Parameters and lengths are in trap units 
[aho = ^yh/{mulho)]■ 



and the velocity v is directed along the x— axis. Further, 
as in the 2D^ case, the harmonic trap is cut off by a ID 
HWBP such that 



ip{~Lx;t) = (piLx;t) = 0. 



(13) 



Vix-,t) = -a;2 + Aexp[-I3{x - ut)^], (12) 



The goal is to provide additional support for our ar- 
guments made about the nature of the observed 2D^ 
SIMW stripes. 

Figure 10 demonstrates the density |(/3(x;i)|2 at a 
time t = 16.1 for TV = 5, A = 10, /3 = 3, and w = 2. 
One can see an SIMW emitted to the left of the cen- 
tral BEC density travelling along the —x direction. The 
sharp dip in the density near x = — 10 is a dark soliton 
excited by the moving obstacle. The DW was able to 
climb up the external potential barrier up to a distance 
of ^ 20 trap lengths from the center of the trap! Amaz- 
ingly, these waves are largely energized, enough to travel 
that far from the center of the trap. Given for example 
a distance of a; = 20a/io, one estimates this wave must 
have an energy of at least 200hiLiho- Hence, a question 
arises whether the moving obstacle is indeed the mech- 
anism which would impart such a high energy to these 
excitations. In fact, if one were to compute the "mass" 
of the obstacle and calculate the associated kinetic en- 
ergy, one would then verify this fact. Thus, let us make 
the following estimation: If the "mass" of an obstacle 
is m and its velocity v, then by assuming the obstacle 
imparts all its kinetic energy to the outgoing ID^ dis- 
persive waves from the trap center, these waves need an 
energy of {l/2)mv'^ ~ 200fiu!fio to climb up the external 
harmonic barrier a distance x = 20aho, before getting 
reflected back from the box potential. Accordingly, if 
V = 2, then m ~ IOOto^, where m^ is the boson mass in- 
corporated in the trap length aho = \/^/{'m'B^ho)- That 
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FIG. 11: Dynamics of the energy (E{t)) [Eq.([8|] of the system 
in Fig. [l] [solid circles, initial condition (i)]. The open circles 
are for Fig. |14| below with initial condition (ii). Energies are 
in trap units (hujho)- 



is, the obstacle posesses a large mass equivalent to the 
mass of ^ 100 bosons, all moving at w = 2. Given the 
latter mass, one can conclude that the moving obstacle is 
able to impart very high energies to the system, exciting 
dispersive waves able to surmount the external potential 
barrier to extremely large distances! 



E. Energy Dynamics 

In this section, we display the time-dependence of the 
energy functional [Eq.Q] for the systems considered in 
the previous sections. The goal is to present evidence for 
a possible existence of solitons inside the 2D^^ trapped 
BEG excited by a moving obstacle, based on the fact that 
the time-dependence of the energy is oscillatory for soli- 
tons [34] • First, we present Fig. 



11 which displays the 
energy dynamics for the system of Fig. IT] [solid circles, 
initial condition (i)]. The solid circles are for Fig. 14 
below with initial condition (ii) added here for the pur- 
pose of comparison. The energy E{t) [Eq.(|8])] fluctuates 
in time about a seemingly stable time average; this is 
indicative of the presence of solitons in the system [53] . 
It can also be argued that the system is emitting and 
reabsorbing energy. 

Figure [12] displays the energy dynamics for the 
anisotropic system of Fig. [3| and Fig. [13] the energy dy- 
namics for the systems in the different trapping geome- 
tries of Figs. 7 9 with pi = P2 = 2.2 (open circles). 
Pi = P2 = 3.0 (solid circles), and pi = P2 — 4: (open tri- 
angles), respectively. Fig. 12 reveals a more rapid and 



spikier behavior in the energy oscillations than the other 
two figures. The system in an anisotropic trap thus tends 
to emit and absorb energy at a faster rate than in an 
isotropic trap. Fig. JTS] displays the effect of the trap 



FIG. 12: As in Fig. [TT] but for the system in Fig. [S] Energies 
are in units of the trap Tiujho- 



curvature on the energy oscillatory patterns. Note that 
the oscillatory behavior becomes more irregular as the 
power of the trap is increased via pi and p2- In fact, the 
energy seems to become almost stable with time when 
t > 4 for pi = p2 = 2.2 and 3. 

Next we proceeded to determine the role of the mov- 
ing obstacle and the nonlinearity of the TDGPE in the 
observed energy oscillations. We therefore checked the 
time-dependence of the energy, Eq.([8) for two cases: 1) 
in the absence of the obstacle and the presence of the 
nonlinearity and 2) vice versa. In case (1) the system 
reached equilibrium immediately after the initialization 
process and stopped evolving. No SIMWs were observed 
whatsoever as the BEG is neither excited nor is the mag- 
netic trap switched off to allow its expansion. The energy 
and density pattern stabilized with time as well. No en- 
ergy oscillations were observed. However, in case (2), 
SIMWs were generated. The system and energy were 
still evolving after the initialization process and the en- 
ergy displayed oscillations. Therefore, the main source 
for the observed oscillations is the moving obstacle. Gon- 
sequently, the moving obstacle not only generates the 
2D^1 SIMW stripes observed in Figs, [ij [s) [e) and [tJIq 
but also solitons inside the central part of the BEG. 



F. Alternative initial conditions 

Here, we demonstrate the effect of using the alter- 
native initial condition (ii) presented in Sec. II E This 
is when the laser potential (obstacle) is initially kept 
stationary, while its intensity is gradually increased up 
to a certain chosen maximum. After this maximum is 
reached, it is set into motion. The resulting density at 
3 is shown in Fig. 14 which is practically the same as 



Fig . [T| b ut using the alternative initial condition (ii) in 
Sec |IIE| The previously-observed SIMW pattern is not 






tq 4 




PFT{kj, ky] t 



M-- 


= 5, 


A 


= 


30 


. ^ = 1 


13 = 


3, 


V = 


= 2 


, t 


= 3 


vi -- 


= P2 




2 


1 


0.0115 

0.011 

0.0105 

0.01 

0.0095 

0.009 

0.0085 

0.008 

0.0075 

0.007 



10 12 14 16 18 20 




ky (aho) 



FIG. 13: Dynamics of the energy {E{t)} [Eq.(|8|] for various 
trapping geometries: open circles: pi = P2 ~ 2.2, system of 
Fig. W\ solid circles: pi — P2 ~ 3, system of Fig. [S] open 
triangles: pi = P2 = 4.0, system of Fig. |9] Energies are in 
units of the trap hojho ■ 
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FIG. 14: As in Fig. [T] but using the initial conditions (ii) in 
SecHlE] 



obtained at i = 3, even at later times (not shown). This 
result holds similary for the anisotropic case tackled ear- 
lier in Fig. [Sj Consequently, one can argue a very pe- 
culiar result, which shows that generation of SIMWs is 
dependent on the initial conditions of an excited BEC. 
We would also like to note that the energy of this sys- 
tem displays an oscillatory behaviour as well (as shown 
in Fig. 11 (open circles)). Henceforth, the observed en- 



ergy oscillations are not a result of specific initialization 
conditions. 



G. Dynamics of momentum density 

Here we try to detect the condensate in order to ex- 
plore its dynamics. For this purpose, the time dependent 
momentum density is computed from the Fourier trans- 



FIG. 15: Momentum density distribution of the system in Fig. 
IT] obtained by the Fourier transform Eq.([9|. Wave vectors k^ 
and ky are in units of the trap (a^J )• 
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FIG. 16: As in Fig. [T] but for the system of Fig. [3] 



form of the spatial density p{x,y;t), as given by Eq.(l9]). 
Figures |15|17| display the momentum density 
PFT{kx,ky;t) of Figs, fll ^ and l7| respectively, and at 
the indicated times t. In all the latter momentum den- 
sity figures, there exists a Gaussian-like, zero-momentum 
density peak centered at {k^^ky) = (0,0). The smooth- 
ness of this peak and the surrounding density in fc— space 
indicate that the BEC excitations in Figs. [T5p7 occur 
in a continuous band of {kx,ky) states. In fact the 
regime around the central BEC peak, so to speak, dis- 
plays a nonzero momentum density surface. By inspect- 
ing Fig. [15] closely, one can observe some rectangular- 
like momentum-density wave structure surrounding the 
central BEC peak. In other words, the effects of the 
geometry-artefact (see Sec. IV B I are manifested even 



in the fc— space of the current system. Fig. 16 shows no 
specific structure around the central peak, but rather a 
smooth momentum density distribution. Nevertheless, 
and upon a close inspection of the momentum density 
along the k^ = ±3 lines in the neighborhood of ky = 0, 
some depression in the density is observed as indicated 
by a slight bending in the plane. Fig. [17] on the other 



10 



PFT{kx, ky] t) 



0.015 
0.014 
0.013 
0.012 
0.011 
0.01 



M-- 


= 5. 


A 


= 30 


, K=l 


P = 


3. 


V = 


= 2. t 


= 2.4 


Pi-- 


= P2 




2.2 

J 


0.016 

0.015 

0.014 

0.013 

0.012 

0.011 

0.01 

0.009 




kx (ct/io ) 



FIG. 17: As in Fig. [T] but for the system of Fig. [7] 



hand displays some valleyed structure in the vicinity of 
the zero-momentum peak. 

Throughout the evolution of the system, the cen- 
tral BEC peak neither disappears nor does it ex- 
pand significantly, although it fluctuates in amplitude. 
Consequently, the zero-momentum condensate density 
PFTiO,0',t) displays oscillatory patterns shown in Fig. 
18 corresponding to the latter systems, 
circles) k 

2; 



4, Pi 



1, Pi 

P2 = 



P2 

Fig. 



2; Fig. 16 



17 



Fig. 15 



(open 



(solid circles) 
(open triangles) k = 1, 



Pi = P2 = 2.2. The osciUations in /Oft(0, 0; t) are a strong 
indication to particles being excited out of the conden- 
sate as pft{0, 0; t) decreases, and particles "falling back" 
into the condensate as pft{0, 0; t) rises again. Therefore 
the presence of SIMW phenomena seems to promote ex- 
citations out of and deexcitations back into the conden- 
sate. For further visual demonstration of the momentum 
dynamics of the system, see the movie in the supplemen- 



tary material section [36], which corresponds to Fig. 15 
Further, one must note that the condensate oscillations 
are substantially damped in the range 8 ~ i ~ 18 for 
Pi=P2-2.2. 



IV. DISCUSSION 

We now set out to find an explanation for the cross- 
like 2D^^ SIMW stripes observed in Figs. ^ [ej and [t] 
In addition, we discuss briefly the work of Ruostekoski et 
al. ISIi, the analogy with ID"' dispersive waves, backflow, 
and a possible existence for vortices and solitons. 



A. Work of Ruostekoski et al. |31j 

Ruostekoski et al. [311 explored the self-interference 
of BEC matter-waves in hard-wall traps of several ge- 
ometries. Their work is very much related to ours and it 
provides significant results. According to these authors, 
the hard walls can be realized by the application of blue- 
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FIG. 18: Dynamics of the zero-momentum density p(0, 0; r) 
for the systems of Fig. [llwith k = I, pi = p2 = 2 (open cir- 
cles); Fig. Is] with K — i, pi = p2 = 2 (solid circles); and Fig. 
It] with K = 4, pi = p2 = 2.2 (open triangles), respectively. 



detuned far-off resonant laser-light sheets. Their goal 
was to study the evolution of a BEC upon the reflection 
of these matter-waves from the hard-walls. Particularly, 
the nonlinear Talbot effect was investigated in a ID hard- 
wall trap, where solitonlike structures were obtained. It 
was noted, that the nonlinear Talbot effect of a coherent 
matter-wave field is analogous to the Fresnel diffraction 
effect in optics. Similarly, we can argue that the cross- 
like matter- wave pattern observed in our Fig. [l] above, is 
also analogous to the Fraunhofer diffraction pattern [37] . 
In general, they initially confined the system in a mag- 
netic field trap added to the hard-wall potential. Upon 
switching off the magnetic trap, the BEC expanded and 
was refiected from the hard-wall boundaries causing a 
self-interference pattern. Ruostekoski et al. conducted 
these investigations in (i) a ID box with an isotropic 

magnetic trap; (ii) a 3D^ highly anisotropic cigar trap 
with hard-walls introduced by laser-light sheets in the 
axial direction of the trap; (iii) a 2D circular boundary, 
and (iv) a 2D square box, both of them with isotropic 
traps. In (i), a complex self-interference pattern is ob- 
tained in the time evolution of the BEC density, charac- 
terized by canals corresponding to the evolution of soli- 
tonlike structures. This is where the Talbot effect was 
observed. The "density holes" observed in the density- 
dynamics corresponded to a gray solitonlike structure. 
A similar solitonlike structure is obtained in the second 
case (ii) of the cigar trap, though with some differences. 
In case (iii), the system is investigated for two subcases: 
initially symmetric, and with an initial broken symmetry. 
In the former ring solitons are observed, in the latter the 
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solitonlike structures break and form vorticity. In (iv), 
the reflections of the matter-waves result also in soliton- 
like structures. 

The main difference between our work and that of Ru- 
ostekoski et al. is that we do not switch off the magnetic 
trap, but rather we excite the BEC by a moving laser- 
potential. In addition, we needed to apply certain initial 
conditions in order to cause an excitation of dispersive 
waves. Another difference is, whereas the latter authors 
observe the Fresnel diffraction effect in the 2D^ square 
box potential, we are able to obtain cross-like matter- 
wave patterns analogous to a Fraunhofer diffraction pat- 
tern. In addition, we would like to emphasize that we 
have focused here on the energy and momentum dynam- 
ics of the BEC matter-wave field. 



B. Artefact of the Geometry 

According to Ref . [5T] . then, the rectangular pattern 
of the waves observed in Figs. [lj[3)[6) and[7]is due to the 
nature of the boundary conditions imposed on the wave 
function at the edges of the square mesh during the sim- 
ulation. Our boundary conditions Eq.([6]) force the wave 
function to vanish there, thus mimicing an impenetrable 
hard wall. According to Eq.([7|), one can see that into this 

box there is embedded the 2D^ harmonic oscillator (or 
PL) trap, both of which share the same center. The hard 
walls of the box are parallel to the coordinate axes x and y 
and the harmonic (or PL) trap ends where the hard walls 
begin and cannot reach beyond. Yet, the motion of the 
obstacle energizes the system to such an extent producing 
waves capable of reaching the edges of the square mesh, 
but not farther than this. At the hard walls the wave 
function of the system vanishes and —possessing enough 
energy— it is reflected from the HWBP. As a result, a self- 
interfering matter-wave pattern forms, as waves outgoing 
from the trap center interfere with "incoming" waves re- 
flected from the hard walls. Inspect also carefully our 
movie in the supplementary section [35j . The perpen- 
dicular 2D^^ SIMWs observed are a result of a trapping 
geometry-artefact arising from the inner curvature of the 
magnetic trap (PL trap), the HWBP, and the finite width 
of the central BEC peak from which the outgoing DWs 
are emitted. The presence of a moving density dip due 
to the obtacle is largely responsible for the presence of 
this cross-like SIMW pattern. Further, one can argue 
similarly to Ref. [SI] that the latter cross-like patterns is 
a result of an initial symmetry breaking of the system 
caused by an initial slight displacement via the moving 
obstacle. On the other hand, the pair of parallel stripes 
observed in Fig. [7] is most likely due to the emission of 
a pair of nonconcentric circular DWs from the central 
BEC peak, each of which is reflected from the hard wall, 
creating its own SIMW pattern. 

Now when the trap is strengthened by increasing its 
power exponents, i.e., by setting pi and p2 larger than 



2.2 in Eq.Q, the expanding BEC waves do not reach the 
HWBP anymore and rather decay by tunneling through 
the edges of the PL trap. As a result, no SIMW patterns 
are formed in Figs. [8] and [9] 



al 

C. Analogy with ID dispersive waves 

We consider an investigation by Hakim [T7], which 
carries some similar features to ours. He explored the 
flow of a nonlinear Schrodinger fluid past an immobile 
obstacle in a one dimensional homogeneous Bose gas and 
found that the obstacle repeatedly emitted gray solitons 
propagating downstream, i.e., opposite to the motion of 
the fluid. At the same time he observed upstream propa- 
gating dispersive waves, i.e., in the direction of fluid flow. 
We can, therefore, explain the initial emission of waves 
from the center of the trap, and before they are reflected 
from the HWBP, based on an analogy with Hakim's in- 
vestigation. Hakim [T7] demonstrated the excitation of 

dispersive waves (DWs) in a ID^ uniform Bose gas in 
the direction of obstacle motion or fluid flow. Basing 
on Hakim's results, we could then argue that the waves 
emitted from the trap center in our 2D system are sim- 
ilarly circular DWs. Consequently, an outgoing circular 
DW interferes with an "incoming" linear DW reflected 
off the box-potential boundaries. However, the cross-like 
pattern observed is a peculiar result which, to the best 
of our knowledge, has not been reported elsewhere. 



D. Vortices and solitons 

On inspecting Fig. [T] closely, one can see two cylindri- 
cal, mantel-like walls surrounding the central BEC peak. 
The region between each "mantel" and the central BEC 
peak displays a deep density depression which is reminis- 
cent of crescent vortex solitons [16J. There is a mantel 
in the front of the central peak (in the region y < 0) and 
in the back (at y > 0). In addition, the pair of "holes" 
which was indicated to earlier in Sec |III A{ is a vortex 
anti-vortex pair similar to the one observed by Jackson 
et al. IJ except that in our pair the vortices are much 
smaller. It is anticipated that vortex excitations could 
play a role in exciting DWs in 2D systems. Indeed, the 
oscillating feature of the energy dynamics shown in Figs. 
|lip3| is a manifestation of solitonic presence inside the 
trap. According to Parker et al. [31], the energy of a soli- 
ton is oscillatory; this supports our previous conclusion 
regarding the possible presence of vortex solitons at some 
instant of time in Fig. [Tl 
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E. Suppression of dispersive-waves and soliton 
generation 

Dispersive waves are suppressed in an anisotropic trap 
along the strongly-confining direction, as displayed in 
Fig. [3] whereas they are still emitted along the major 
axis of the elongated trap. The dispersive waves do not 
have enough kinetic energy to climb up the external po- 
tential barrier along the strongly-confining direction (the 
minor axis of the trap) . 

Next to this, upon increasing the curvature of the 
isotropic confinement, as in Figs. [7][9J the number of 
solitons inside the trap increases significantly. The sup- 
pression of DWs prevented the energy necessary for soli- 
tonic excitations to be carried away by dispersion and 
be expended on climbing up the external potential bar- 
rier. One can conclude that a larger number of solitons 
can be generated in a strongly-confining PL trap by a 
moving obstacle. 



F. Backflow 

According to Feynman and Cohen (1953) (FC) [55], 
the motion of an obstacle in a fluid causes backflow of 
the atoms in the neighborhood of the moving obstacle. 
Referring back to the FC theory, an excited atom i in a 
strongly interacting Bose fluid such as liquid ''He is sur- 
rounded by adjacent atoms j that correlate with atom i 
and are thus able to flow around it. The correlation func- 
tion can be found in Ref. [SF . As such, one could further 
note that collective excitations in our current system gen- 
erate waves travelling away from the center of the trap 
which are later reflected from the 2D^ box-potential. We 



further believe that in the initial condition (i) in Sec II E 



the gradually introduced laser potential induces a back- 
flow pattern which, at the end of the initialization stage, 
interacts with the backflow pattern of the secondary laser 
potential introduced in the transient stage of the simu- 
lation. This is in a manner consistent with the backflow 
physics discussed in an earlier investigation by Ghassib 
and Chatterjee [39] . where they computed the interaction 
potential between the backflow patters resulting from the 
motion of two impurities in a non- viscous fluid in various 
geometries. 



G. Future w^ork 

The numerical analysis of the current system con- 
ducted by applying the Crank Nicolson method to the 
2D TDGPE lead us to speciflc information on the dy- 
namics of this system as elaborated on in Sec. |III| How- 
ever, other information regarding the quantum hydro- 
dynamics of this system is of importance and needs to 
be investigated. Particularly the phase of the wavefunc- 
tion could reveal further information about the coherence 



properties. For this purpose, the Madelung transforma- 
tion (MT) [23] can be utilized in obtaining two analytic 
equations arising from the real and complex part of the 
MT. In this regard, the dynamic phase of the wavefunc- 
tion reveals information about the dynamic superfluid 
velocity field in the presence of the moving obstacle po- 
tential. In the future, we will investigate also trapped 
BECs excited by an attractive (red-detuned) laser pote- 
nial. 



V. CONCLUSION 

In summary, then, we have demonstrated the cre- 
ation of highly energetic, two-dimensional, self interfering 
matter-wave (SIMW) stripes inside a 2D^ BEC confined 
by a power-law (PL) potential trap cut off by hard-wall 
box-potential (HWPB) boundaries. The split-step Crank 
Nicolson method was used to numerically solve the time- 
dependent Gross-Pitaevskii equation (TDGPE) using a 
code previously developed by Muruganandam and Ad- 
hikari [30]. It was found that the SIMWs are generated 
by the interference of circular dispersive waves (DWs) 
emitted from the center of the trap and their refiected 
wave fronts off the box-potential walls. Four mutually 
perpendicular SIMW stripes are observed, two of them 
parallel to the obstacle's motion along the j/— axis and 
the others perpendicular to y. As we have indicated, the 
DWs are chiefly excited by an obstacle moving inside the 
2D BEC under speciflc conditions. These conditions in- 
clude an interaction strength of the order of A/" ~ 1 and 
harmonic confinement, or a PL trap, given by Eq.Q with 

Pi,P2^2.2. 

The suppression of dispersive waves by strong confine- 
ment results in a buildup of energy inside the BEC, which 
eventually leads to a larger number of soliton excitations 
than when this energy is expended on DW excitations. 
By Fourier transforming the spacial density distribution, 
the momentum density dynamics displayed a nonvanish- 
ing, oscillating zero-momentum condensate density. As 
the strength of confinement was increased, the latter os- 
cillations began to disappear in a certain time range, as 
revealed in Fig. [Tsjfor pi = P2 = 2.2, indicating that the 
condensate fraction began to stabilize somewhat in that 
regime. Particularly, the condensate density oscillations 
signal the excitation and deexcitation of particles from 
and to the zero-momentum state. Further, a vortex anti- 
vortex pair was obtained in the vicinity of the moving 
laser obstacle similar to, but smaller in size, to those of 
Jackson et al. [Tj. 

In essence, then, we presented in this regard a phe- 
nomenological examination of these systems in analogy 
to earlier examinations by Hakim 17J and He et al. [16j . 
The excited waves were classified as dispersive based on 
the earlier ID^ investigation of Hakim [17], who found 
that a nonlinear Schrodinger fiuid fiowing past a station- 
ary obstacle can emit dispersive waves. By considering 
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the inverse of this argument, i.e., a moving obstacle inside 
a stationary BEC, we arrived at the latter classification. 
In addition, following somewhat Hakim's investigation, 
we conducted ID^ simulations for a BEC trapped in a 
harmonic potential and were able to observe density os- 
cillations excited in the BEC. These oscillations shown 
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